Spectral properties of incommensurate charge-density wave systems 
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The concept of frustrated phase separation is applied to investigate its consequences for the 
electronic structure of the high T c cuprates. The resulting incommensurate charge density wave 
(CDW) scattering is most effective in creating local gaps in k-space when the scattering vector 
connects states with equal energy. Starting from an open Fermi surface we find that the resulting 
CDW is oriented along the (10)- and (or) (Ol)-direction which allows for a purely one-dimensional 
or a two-dimensional "eggbox type" charge modulation. In both cases the van Hove singularities 
are substantially enhanced, and the spectral weight of Fermi surface states near the M-points, tends 
to be suppressed. Remarkably, a leading edge gap arises near these points, which, in the eggbox 
case, leaves finite arcs of the Fermi surface gapless. We discuss our results with repect to possible 
consequences for photoemission experiments. 



I. INTRODUCTION 

Striped phases are an important issue in the dis- 
cussion on physical properties of the high-T c materials 
|0,|). Incommensurate magnetic peaks displaced from 
the antifcrromagnetic wave-vector along k x and k y have 
been observed by inelastic neutron scattering in both 
La^Sr.CuCU @-§ and YBa 2 Cu 3 7 _ ;c (YBCO) [§, 
whereas in the Bi 2 Sr 2 CaCu208 material the incommen- 
surate magnetic peak position has not yet been resolved 
M. However, incommensurate magnetic scattering may 
result from either some kind of spiral phase (involving a 
modulation of the transversal spin components only) or 
from longitudinal spin fluctuations where in this case a 
strong coupling to the charge is expected. Evidence for 
cooperative charge- and spin scattering is based on the 
observation of domain walls in Laig-^Ndo^Sr^CuO^ |8J] 
where the low temperature tetragonal lattice structure 
and the filling close to 1/8 are suited to pin the den- 
sity fluctuations, giving rise to a static CDW phase as 
revealed by commensurate ionic shifts observed in neu- 
tron scattering. From the temperature dependence of 
the charge- and spin-order peaks one can further con- 
clude that the stripe order is driven by the charge rather 
than by the spins. Additional support for charge induced 
stripe correlations is provided by recent neutron scatter- 
ing studies of the phonon dispersion in YBCO M . These 
measurements show a large broadening in the spectrum 
at wave vectors which are twice the magnetic incommen- 
surability. Moreover, the broadening in the phonon spec- 
trum occurs at higher temperatures as the incommensu- 
rable magnetic peaks indicating again that charge fluc- 
tuations are the 'driving force' behind incommensurable 
magnetic scattering 



The formation of stripe correlations in these com- 
pounds can be well understood within the concept of 
frustrated phase separation [pi where a phase separation 
instability is prevented by long-range Coulomb interac- 
tions [j[3|. As a result the long- wavelength density fluc- 
tuations associated with phase separation are suppressed 
in favor of shorter-wavelength density fluctuations, giv- 
ing rise either to dynamical slow density modes flCf| o r 
to incommensurate charge density waves (CDW) |p.|,[Ll|- 
Concerning the pairing mechanism the theoretical pro- 
posals differ with regard to the role of the stripe corre- 
lations on the superconducting pair formation. Whereas 
some theories consider the formation of quasi 1-d 'rivers 
of charge' as a necessary environment for the quasiparti- 
cles to pair or to enhance T c [Q,|l2[ others take the fluc- 
tuations associated with the stripe instability itself as a 
pairing mechanism [GJ . 

The proximity to a stripe instability tuned by doping 
and temperature allows to interpret many properties of 
the cuprates as due to a Quantum Critical Point (QCP) 
JUJ located near optimal doping HlJJf-fOJ. Within this 
scenario the singular scattering induced by the critical 
fluctuations would be responsible for both the anoma- 
lous normal-state properties and the large superconduct- 
ing critical temperatures. Then the phase diagram of 
the cuprates is partitioned in the (nearly) ordered, the 
quantum critical, and the quantum disordered regions 
corresponding to the under-, optimally, and over-doped 
regions of the phase diagram of the cuprates. 

The relevance of the stripe phase and of the stripe fluc- 
tuations should find an experimental confirmation (or 
a disproof) in the analysis of the single-particle spec- 
tral properties. Taking the point of view that incom- 
mensurate CDW scattering contributes to determine the 



electronic properties of the high-T c materials the ques- 
tion arises how the electronic structure around the Fermi 
energy is affected by this scattering and if this can be 
detected by angle-resolved photoemission spectroscopy 
(ARPES) pa]. In fact, there are several features which 
seem to be common to all p-type cuprates (see e.g. Ref. 
jl9| and references therein). In particular, the disper- 
sion displays an extended saddle-point at (ir/a, 0) (here 
and in the following we will consider unitary lattice spac- 
ing a = 1) with an energy close to Ep, which has been 
observed in Bi2212 Q as well as in YBaCuO Q sam- 
ples. Furthermore, the Fermi surface of Bi2212 is very 
strongly nested with a nesting vector of approximately 
0.9 (it/ a, it /a). The resulting picture for the Fermi sur- 
face therefore consists of tube-like structures around the 
M-points which are connected by more or less straight 
pieces. It has been shown in Ref. [g^] that many of these 
unusual Fermi surface features can be explained by as- 
suming a disordered stripe phase to be present in the 
high-T c copper oxides. The fact that the charge carriers 
are exposed to a stripe potential mirrors in the Fermi sur- 
face through one-dimensional characteristics which pre- 
dominantly appear around the M-points. 

A further fascinating Fermi surface feature of the 
high T c superconductors, which has been detected in 
ARPES measurements concerns the pseudogap phe- 
nomenon which is observed in optimally and underdoped 
bismuth compounds below a doping dependent tempera- 
ture T* > T c |§-|2e|. In underdoped samples of Bi2212 
p3| , p4| there is a qualitative change in the electronic 
structure with respect to the optimally doped system al- 
ready above T c , with the appearance of so-called leading- 
edge gaps around the M points. In these k-space regions 
there is no indication of a quasiparticle peak, and the 
spectra are dramatically broadened and shifted to higher 
binding energies with a maximum at 100-200 meV. As a 
consequence large portions of the Fermi surface around 
these points are not visible whereas along the r — ► X 
directions the spectral lincshapes are similar to those of 
the optimally doped samples. Similar results where ob- 
tained in the one-plane material Bi2201 |2^] showing that 
bilayers are not essential to the pseudogap phenomena. 

A common explanation for the features described 
above is based on the existence of pre-formed pairs al- 
ready above T c , since a continuous evolution of the spec- 
tra from the normal to the superconducting state |2q ] 
was apparently observed. In particular, a superconduct- 
ing d-wave gap is in accordance with the symmetry of 
the pseudogap. Moreover, for materials with kp^o < 10 
the existence of two temperature scales for pairing and 
phase coherence is a rather natural consequence. How- 
ever, this simple scenario is put in jeopardy by various 
experimental findings. First of all the rapid increase of 
T* by decreasing doping is to be contrasted with the 
much more moderate increase of the maximum gap at 
low temperature. (T* doubles when T c goes from 90 K 



to about 75 K in Bi-2212, whereas in the same doping 
range, A(T = 0) around the M-points changes by thirty 
per cent at most). This suggests that a substantial tem- 
perature dependence of the pairing potential is present 
in this doping range. Furthermore recent ARPES exper- 
iments ]^9| show that the pseudogap above T c gradually 
appears at the M points leaving extended segments of 
the Fermi surface gapless. These segments continuously 
shrink to a point-like node only below the superconduct- 
ing critical teperature T c . Needless to say that these 
behaviors do not correspond neither to the usual BCS 
gap opening process nor to the phase locking of uniform 
preformed pairs. Instead it might find an interpretation 
along the line of coexisting pairs and fermionic quasipar- 
ticles H]|l|]. 

In front of this unsettled issue, one may ask whether 
the rich and unexplained behavior of the pseudogap phe- 
nomena could result from an electron-hole pair scattering 
above T c rather than from pair correlations alone. This 
line of argumentation has been followed in several anal- 
yses based on the scattering between fermions and col- 
lective spin excitations ]32|-p4[. A quite similar outcome 
would instead arise from the quasiparticles being scat- 
tered by incommensurate charge fluctuations p7[ . The 
fact that a substantial part of the spectral features of 
the optimally and underdoped cuprates can be explained 
in terms of particle-hole scattering is also supported by 
recent experiments |2q] measuring the Fermi surface by 
sequential angle-scanning photoemission. These results 
clearly show the missing segments of the Fermi surface 
around the M-points with a shape, which can be inter- 
preted in terms of scattering of the quasiparticles with 
quasicritical mixed spin and charge fluctuations pj| . 

The idea that the pseudogap could (not only) arise 
from pairing in the particle-particle channel, but also 
from different scattering mechanims (like CDW fluctu- 
ations) finds an additional support in the recent analysis 
of Ref. |3£| , where the experimental results on the low- 
temperature penetration depth suggest that the doping 
behavior of the low-temperature gap near the nodes in 
the (1,1) direction is similar to what expected from stan- 
dard d-wave BCS theory with the gap scaling with T c . 
This is in contrast with the behavior of the pseudogap 
measured at the k-points around (±7r/a, 0) and (0, ±7r/a) 
in ARPES experiments below T*, which increases with 
decreasing T c . This suggests the possibility that at high 
temperature (T* > T > T c ) the gap is due to CDW 
fluctuations scattering in the particle-hole channel. In 
this regard it was shown in Refs. [jjllllj that the incom- 
mensurate CDW scattering resulting from the competi- 
tion between phase separation and long range repulsive 
Coulomb forces, occurs near some scattering vector q c 
not related to the Fermi vector k^. It was also shown 
within a standard large-N approach |l],pT[ that the ef- 
fective scattering amplitude close to q c in the particle- 
hole channel is strongly attractive and similarly struc- 



tured in momentum space in the particle-particle channel 
also, thus resulting in substantial pairing in the d-wave 
channel. Within this scenario, the strong momentum- 
dependent effective interaction is responsible both for 
(preformed) pairing and for pseudogap formation due to 
dynamical charge modulation. Therefore searching for 
spectral features arising from scattering mechanisms act- 
ing in the particle-hole channel is not necessarily in con- 
trast with the interpretation of other features in terms of 
pairing in the particle-particle channel. 

In the present paper we are interested in the question 
if incommensurate CDW scattering, resulting from the 
competition between phase separation and long range 
repulsive Coulomb forces, can account for some of the 
above anomalous features observed in photoemission ex- 
periments. This problem was partly tackled in Ref. |27[| 
within a perturbative treatment of quasicritical spin and 
charge fluctuations. This perturbative approach is unfor- 
tunately not feasible in the underdoped region, deep in- 
side the stripe phase. Therefore we perform here a mean- 
field analysis, which should capture the non-perturbative 
character of well-formed stripes. Of course, the dynam- 
ical character of the stripe fluctuations will be missed. 
Nevertheless, at short distances and time as those probed 



in photoemission experiments, this description should be 
appropriate. 

We emphasize once again that the complexity of the 
physics of the cuprates, particularly of the underdoped 
ones, involves spin and Cooper as well as charge density 
fluctuations. Our aim is therefore not to account for all 
the spectral properties of these materials, which likely 
arise from the interplay of all these mechanisms. We 
rather conservatively aim to point out that some of the 
prominent (and puzzling) features of the cuprates could 
find a natural interpretation by assuming charge quasi- 
ordering. 

The paper is structured as follows. In Section II the 
model and the general formalism are presented. Section 
III contains the physical results for the Fermi surface and 
the spectral densities, which are discussed and summa- 
rized in the conclusive Section IV. 



II. FORMALISM 

Starting point is the Hubbard- Holstein model with an 
additional in-plane long-range Coulomb interaction 
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where Ci a {c ity ) destroys (creates) an electron with spin 
a at site i and Ai{A[) destroys (creates) a local Holstcin- 
type phonon at site i. The summation over nearest- and 
next-nearest neighbour sites is indicated by < i, j > and 
<< i,j >>, respectively. 2 CT n, CT = J2<j c L c i<? is tne lo_ 
cal electron density and its Fourier transform is given 

by Pq = Y,ka c \+qa c k<y The l ast term in Eq.(@) de- 
scribes the Coulomb interaction between electrons on a 
two-dimensional, square lattice (lattice constant a in x- 
and y-direction) , which is considered as a plane embed- 
ded in the 3-dimensional lattice (plane distance d in the z- 
direction). The dielectric constants in the plane and per- 
pendicular to it are given by eii and ej_, respectively. The 

Coulombic coupling constant is V c = 2£ e ai . On the z=0 
plane the momentum dependence of the Coulomb poten- 
tial is found to be G(q x ,q y ) — [e||/e^(a/d) 2 ][cos(ag x ) + 
cos(aq y ) — 2] — 1. As usual, the sum in the Coulombic po- 
tential does not include the zero-momentum component, 
since we are supposing that the diverging q — interac- 
tion between electrons is canceled by the contribution of 
a uniform positively charged ionic background. 

Since we are interested in the limit of strong local re- 
pulsion, we take the limit U — > oo, which gives rise to the 
local constraint of zero double occupancy ^2,a ni < J — 1- 



To implement this constraint we use a standard slave- 
boson technique |39,Ej by performing the usual substi- 



tution c]„ — > c]„bi, Cia — ► Ciab] and implement the lo- 



cal constraint by a local Lagrange multiplier A^. The 
model can first be solved in the mean field approxima- 
tion by setting the b\ and Ai bosons to their constant 
self-consistent values bo and Ao, respectively. At this 
stage the phonons and the Coulombic interaction de- 
couple from the fermionic quasiparticles for an homo- 
geneous mean-field solution. The system then results 
in free quasiparticles with a shifted chemical potential 
P = Po ~ Ao and a dispersion E k = — 2£& ! e/ £ , where e k = 
[cos(ak x ) + cos(ak y ]+t' /t[cos(ak x + ak y ) + cos(ak x — ak y )] 
and Oq — 5 corresponds to the concentration of doped 
holes. The mean- field value for Ao is determined by 
Ao = AO + {t'/t)Xl = 2tJ2 k f(E k )e k where f(E) is the 
Fermi function. 

The effective interaction leading to scattering between 
quasiparticles arises from the exchange of the bosonic 
fields beyond the mean-field approximation. Working 
within the radial gauge |4fj] one can define a three- 
component field a^ = (5r, SX, <5</>) where S(j) is the lattice 
displacement field, and 5r, 5X are the fluctuating parts 
of the modulus of the b\ -bosons and of the Lagrange 



multiplier A^ respectively. The static bare scattering am- 
plitude in the particle-hole channel can be written as flO] 

r (fc, fc'; q) = -J2 AM ( fc '> -q)D% v {qW{k, q) (2) 



where the vertices A M coupling the fermionic quasiparti- 
cles to the bosons are defined as 

A r (k,q) = -2tbl(e k+q/2 + £k- q /%) 
A\(k,q) = i 
A cf) (k,q) = -2g 

and D^ v (q) denotes the matrix of the bosonic Greens 
function JT1J with all elements being zero except for 

DY = (b 2 /2) [X° Q [sm 2 (q x /2)+sm 2 (q y /2)} 

+ {t'/t)Xl[sin 2 {{q x - q y )/2) + sin 2 (( fe + q y )/2)] ] , 
D r ' x =D^=ibl/2, 
Dp* = ^. (3) 

We finally obtain a Hamiltonian, describing the effec- 
tive interaction between quasiparticles 



ka q 



where the static effective interaction is given by 



(4) 



V q = T(k F , k' F - q :U ; = 0) = U + j q ~^- 
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U = -4E F /5 and 7g = ^[sin 2 (g x a/2) + sm 2 (q y a/2)} + 

^f[sin 2 (( fe - q y )a/2) + sin 2 (( fe + q y )a/2)]. Within a 
Fermi-liquid scheme, V q represents the effective residual 
interaction between the quasiparticles once the screen- 
ing due to intraband particle-hole quasiparticle bubbles 
is disregarded. 

The Hamiltonian (Q) is the starting point of all fur- 
ther investigations. Its interaction consists of a part, 
U + 7g, increasing with |q|, which describes the residual 
scattering of the quasiparticles through the slave-bosons 
exchange. The superposition of the long-range Coulom- 
bic part causes the interaction to exhibit a minimum as 
a function of q. The electron-phonon coupling 2g 2 /loq 
rigidly shifts the potential to lower values. This will in- 
duce an instability in the density-density response func- 
tion at a critical value of q = q c when 1 + V^ c II(q c ) = 0, 
where n(q) is the particle-hole fermionic bubble. This 
instability marks therefore a transition to an incommen- 
surate CDW with scattering vector q c Hj. 

To proceed we will approximate the effective model 
Eq. (H) by a Hartree factorization and assume that the 
symmetry of the system is broken with respect to a given 



vector q c and its multiples. Notice that retardation ef- 
fects arising from the frequency dependence of the total 
scattering amplitude have been neglected at this level. 
This mean-field treatment of Eq.[| finds some justifica- 
tion in investigating static or low-energy properties of 
the quasiparticles close to the Fermi level, whereas this 
approximation is crude in the analysis of dynamical prop- 
erties at frequencies, which are comparable with the typ- 
ical energies of the exchanged boson. 

In order to reproduce the 'unperturbed' Fermi surface 
(i.e. of the overdoped system) of the Bi22I2 compounds 
at optimal doping we take the hopping parameter t' to 
be t' = — 0.45i, which leads to an open Fermi surface 
centered around (±7r/a, ±7r/a). It turns out that within 
a linear response approach as in Ref. jnj, for this kind 
of bandstructure the CDW instability occurs first in the 
(1,0) or (and) (0,1) directions. In fact, a negative value 
of t' enhances the density of states around the M-points, 
which mirrors in an enhancement of the fermionic bub- 
ble II(q) in the (1,0) and (0,1) directions. Thus in a 
RPA-like treatment the divergence of the density-density 
response function preferably takes place for q = q c in 
the (1,0) or in the (0,1) directions. Within a linear- 
response analysis, the instability is forced to occur at a 
given wavevector and the above arguments are in favor 
of the instability taking place with a modulation alter- 
natively in the x or y direction described by the order 
parameter < p q > 1D = J^n < Pg > & qnq */v However, 
a non-linear instability (as the one here described by a 
standard set of Hartree mean-field equations) can also 
occur simultaneously along both directions leading to a 
two-dimensional eggbox type modulation of the charge. 
In this case, the density expectation value is given by 

< p q > £ 99= Y,n < Pi > l S q,riq S + ^q,nq v c }- 

Although we are now dealing with an effective one- 
particle Hamiltonian, it is clear that the system cannot 
be diagonalized for general incommensurate q c vectors. 
We therefore represent q c as 



7r n x rij 

q c = - — ; — 

a m x m 



(6) 



which allows us to simulate the 'incommensurability' 
by increasing the values for both m and rm, respec- 
tively. Once the the scattering term in Eq. ^ is de- 
coupled a la Hartree, for each k-point in the reduced 
BZ, the bare Bloch functions v I / k(«) = ^(k + nq c ) with 
k-vector k, k + q c ,k + 2q c , ....,k + n max q c , are mixed 
by the interaction. n max depends on q c via the con- 
dition k + (n max + l)q c = k. For each k-point, the 
Hamiltonian can be put in a matrix form, which can 
be diagonalized thus finding the linear transformation 
from the bare Bloch states tykin) to the new eigenstates 
^k(n-) = Xim a k( Ti ' m )^k('n)- The system can be di- 
agonalized numerically and the CDW order parameters 
Xn = Vnq a < Pnq c > are determined self consistently. 



Once self-consistency is reached, the eigenvalues Ek(m) 
give the band structure of the system. 

The spectral function in the full BZ then consists of a 
ensemble of weighted delta-functions and is given by 



^k+nq c M = ^al{n,m)S(uj - E^{m)) 



(7) 



It should be noticed that this approach mantains its 
full validity as long as the initial constraint of no double 
occupancy is satisfied. Therefore we will restrict in the 
following to the cases where the CDW order parameter 
is small enough that the local density never exceeds one, 



III. RESULTS 



In this section we will discuss the bandstructure, Fermi 
surface and photocmission spectra resulting from cither 
1-d or eggbox type CDW scattering. For simplicity, in 
the following a planar cell of unitary length (a = 1) will 
be used. Since in the sytems with two Cu02 layers per 
unit cell (as in Bi2212) different stripe orientations in 
different planes or also in different regions of the same 
plane may be realized, we consider also the case of su- 
perimposed (10)- and (Ol)-stripes. On the other hand 
the eggbox-solution, where the symmetry is broken in 
the x and y directions at the same time, can be regarded 
as a model for 1-d stripes fastly fluctuating between the 
(10)- and (Ol)-direction, while maintaining a fixed |q c |. 
Furthermore, having in mind the possibility of dynamical 
CDW scattering we also investigate averaged disordered 
realizations of an eggox charge modulation, where fluc- 
tuations in the absolute length of the scattering vector 
will be considered. 



A. One-dimensional CDW phase 



The results in this section supplement considerations 
in Ref. [p2[ where some effects of 1-d CDW scattering on 
the spectral properties were investigated. Here we ana- 
lyze the consequences of CDW scattering on the band- 
structure and, differently from Ref. |22] , we start from an 
open Fermi surface (induced by a next-nearest neighbor 
hopping term t' as discussed in the previous section). 




FIG. 1. Fermi surface for a 1-d charge modulation near the 
CDW instability where only the first harmonic xi — 0.04 is 
different from zero. Parameters: 8 = 0.2, a — —0.45, u>o = 
0meV,g = 0A6eV,V c = 1.65eV. CDW modulation \q c \ = 
0.79 in direction of V — M. The plot is for temperature 
T — 100A" and the energy window around Ef has choosen to 
be 50meV. Intensities: I > 50%: full points, 10% < I < 50%: 
circles, 1% < I < 10%: small dots. 

Fig. 1 displays the calculated Fermi surface (FS) of a 
one-dimensional (1,0)-CDW for a parameter set close to 
the instability and doping S — 0.2. Parameters such that 
q c pa (jr/Aa, 0) were chosen, corresponding to a charge 
modulation with a period of 8 planar unit cells. Since we 
diagonalize the system in a reduced Brillouin zone (BZ) 
defined by the choice of q c in eq. (||) the defolding of the 
bands leads to a redistribution of spectral weight for the 
FS states. The photoemission intensity for each k-point 
is given by the following integration over an energy win- 
dow around the Fermi energy Ik = J E F _ e dujf(u))A\ i (u>) 
where f(oj) denotes the Fermi function and A-^(uj) is 
the spectral function defined in eq. (R). For simplic- 
ity the intensity range is divided into three ranges for 
high (Ik > 0.5, full circles), intermediate (0.1 < Ik < 0.5, 
open circles) and low (0.01 < h < 0.1, dots) intensity. 

Similar to the results of Ref. |2^] there appear dis- 
placed shadow bands with low intensity which are due 
to q c -scattering processes of the original FS states. How- 
ever, since we started with an open Fermi surface (due to 
a quasi one-dimensional bandshape around (tt, 0), (0, 7r)) 
the k-space regions around the M-points are strongly af- 
fected by the scattering. Moreover, there appears an 
asymmetry between these two areas since around the M- 
point the scattering is along the FS branches whereas at 
M it connects the branches centered around X and Y. 
As a consequence more k-states are reduced in intensity 
around (n, 0) since the number of states with approxi- 



mately the same energy that can be connected by q c is 
larger than around (0, 7r). 
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FIG. 2. (a-c): Cuts of the bandstructure in the full 
Brillouin zone along various symmetry directions for (10)- 
scattering, corresponding to the Fermi surface in Fig. 1. (d): 
Superposition of the M — > X and the M — * X cut, corre- 
sponding to the Fermi surface in Fig. 1. Intensities: / > 50%: 
full points, 10% < I < 50%: circles, 1% < I < 10%: small 
dots. 

To analyze the electronic structure in more detail we 
plot in Fig. 2 the bands with weight larger than 1% in 
the interesting symmetry directions F — > X and M — ► X. 
The intensity range is divided in the same way as in 
Fig. 1 and the same symbols are used. Moving along 
r — ► X (Fig. 2a) one observes a rich gap structure around 
r which diminishes upon approaching the FS crossing 
where besides the main band only the two weak shadow 
bands survive. The multiband features around T are due 
to the fact that at the bottom of the band its slope is 
small and therefore many states with similar energy can 
be connected by the CDW vector. Taking the cut along 



M — > X (Fig. 2b), which is orthogonal to q c , one ob- 
serves that the CDW has induced the formation of a sec- 
ond band which upon approaching X rapidly loses in in- 
tensity. On the other hand the band along the M — > X 
direction (Fig. 2c) displays a CDW gap at the FS cross- 
ing since this direction is parallel to q c and the size of 
the scattering vector is of the same order than the FS 
branch separation at (71", 0). 
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FIG. 3. Photoemission spectra corresponding to the cuts 
in Fig. 3. The fermi function has been added as a back- 
ground. The broadening of the 5-functions is 25meV and the 
temperature T=100K 

In Fig. 3 we show the energy distribution curves corre- 
sponding to the bandstructure cuts of Fig. 2. Along the 
r — > X direction (Fig. 3a) the gapped structure around 
r appears as a broad multipeaked feature which evolves 
into a single peak upon moving towards the X-point. 

For the cut along M — > X one observes two peaks 
crossing the FS at (ir, 0.39) and (n, 0.52) corresponding 
to the two bands in Fig. 2b. The CDW along M -> X 
displays in a peak hrst moving towards the Fermi energy 
and then bending back to lower energies (Fig. 3c). How- 



ever, this gap no longer can be detected when we consider 
a system with superimposed (10)- and (01) CDW scatter- 
ing. The corresponding Fermi surface, bandstructure and 
energy distribution curves are plotted in Figs. 4, 2d, 3d re- 
spectively. Although there is still reduced weight of the 
k-states around (0,7r), (n, 0), which are now equivalent, 
the CDW gap at M is no longer visible in Fig. 3d since 
the averaged bandstructure is dominated by the states 
around the M-point. 
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-3.2 -2.4 -1.6 -0.8 0.0 0.8 1.6 2.4 3.2 
FIG. 4. Fermi surface for the superposition of the (10)- 
stripe of Fig. 1 with the equivalent (Ol)-charge modula- 
tion. Parameters: S — 0.2, a = —0.45, ojo = 0.06eF, g = 
0.46eV,V c = 1.65eV. CDW modulation \q c \ = 0.79. The 
plot is for temperature T = WOK and the energy window 
around E_f has choosen to be 50meV. Intensities: I > 50%: 
full points, 10% < I < 50%: circles, 1% < I < 10%: small 
dots. 

It is also of interest to investigate the band disper- 
sions along the r — > M and r — » M directions. Along 
these directions a substantial enhancement of the van 
Hove singularities have been observed by ARPES experi- 
ments |20||2l]Jlq| and the obvious question arises whether 
this feature can arise from CDW scattering. Fig. 5 re- 
ports in the (a) and (b) sectors the results of our calcu- 
lation starting from an unperturbed bandstructure (indi- 
cated by the diamonds in Fig. 5a,b). It is apparent that 
multiple non-dispersive shadow bands appear along the 
r — > M direction as a result of the CDW scattering, with 
a substantial redistribution of the spectral weight both 
at higher and lower energies. This latter occurrence sug- 
gests that indeed the broadening of quasiparticle spectra 
towards lower energies could experimentally be seen as 
a broadening in fc-space of the area where the band is 
non-dispersive. For the presently considered (10) one- 
dimensional scattering, the same does not hold in the 
r — » M direction, where the shadow bands are strongly 



dispersive and much less numberous (Fig. 5b). As a con- 
sequence the spectra around M would hardly look like 
an enhanced van Hove singularity at low energy. 
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FIG. 5. (a-c) Bandstructure for one-dimensional CDW 
scattering along the (lO)-direction for the same parameters 
as in Fig. 1 of the paper, (a) T — > M — (n,0), (b) 
r — > M — (0, 7r), (c) Superposition of (a) and (b). The 'un- 
perturbed' band is indicated by diamonds, (d) The same cut 
for 2-d eggbox scattering. Here the 'unperturbed' band is in- 
dicated by the full line. Intensities: / > 50%: full points, 
10% < I < 50%: circles, 1% < I < 10%: small dots. 

However, the possibility still remains open that the 
spectra of the real materials can result from the (nearly 
static) superpositions of the spectra in Fig. 5a, b as a 
consequence of the different orientation of the stripes on 
different Cu02 planes or on different regions of the same 
plane. This is depicted in Fig. 5c. Alternatively, if the 
stripes are fastly fluctuating, their effect on the spectra 



could better be mimicked by our (static) eggbox solution 
described below. 



B. Ordered eggbox phase 

In this section we are concerned with a two- 
dimensional charge modulation resembling the shape of 
an eggbox. 

Fig. 6 displays the calculated FS of a regular 2-d CDW 
for a parameter set close to the instability and doping 
5 = 0.2. Similar to Fig. 1 the scattering mostly affects 
the states around (0,7r) and (n, 0), but now of course 
the features at M and M are symmetric. Moreover, the 
region of reduced intensity is much larger than for the 
superimposed (01)- and (10)-CDW in Fig. 4 (note that 
in both plots the strength of the order parameter is the 
same) . 
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FIG. 6. Fermi surface for an eggbox type charge modula- 
tion near the CDW instability where only the first harmonic 
Xi = 0.04 is different from zero. Parameters: S — 0.2, a = 
-0.45, ujo = 0MeV,g = 0A6eV,V c = 1.65eV. CDW modu- 
lation |g c | = 0.79. The plot is for temperature T = WOK 
and the energy window around Ef has choosen to be 50meV. 
Intensities: I > 50%: full points, 10% < I < 50%: circles, 
1% < I < 10%: small dots. 



(less than 10% of intensity) features at (1.6,1.6). It is 
not unconceivable that these multiple crossings of the 
(shadow) Fermi surface(s) have been observed, although 
the common interpretation refers to magnetic scatter- 
ing 1 33 1 . This finding suggests that, together with the 
magnetic scattering, charge scattering could cooperate to 
bring spectral weight at the Fermi surface in a "pocket- 
like" shape §7|. 
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FIG. 7. Bandstructure in the full Brillouin zone corre- 
sponding to the Fermi surface in Fig. 1. Intensities: / > 50%: 
full points, 10% < I < 50%: circles, 1% < I < 10%: small 
dots.(a) r < k < X, (b) M < k < X. 



An additional quite interesting feature is the appear- 
ance of shadow Fermi surfaces in the (1, ±1) directions, 
giving rise to "pockets" . This feature was already hinted 
to in the Fermi surface of Fig. 4, where it was, however 
less pronounced. In particular, here moving from the 
r point towards the X (or the other equivalent points) 
one first finds some spectral weight at the Fermi energy 
for points located around (0.8, 0.8) in fc-space. Moving 
further towards X one meets the main Fermi surface 
with highest intensity around (1.2,1.2), and then weak 



The effect of the 2-d scattering becomes more trans- 
parent in the bandstructure cuts plotted in Fig. 7. Scan- 
ning along the diagonal direction gives a strong splitting 
of states around the T point similar to our previous find- 
ings for 1-d scattering. Moving further towards the X- 
point most of these states rapidly lose in intensity and at 
the FS crossing there is only one main band left together 
with two shadow bands with low intensity in accordance 
with the FS shown in Fig. 6. 
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approaching the FS crossing, due to the rapid loss of in- 
tensity in the scattered bands, there appears a evolution 
into a single peak structure which crosses the Fermi level 
at k = (1.18, 1.18). In the right panel of Fig. 8a we plot 
a comparison of the CDW spectra with the peak evo- 
lution of the homogeneous system for selected k-points. 
From this plot it becomes clear that along the diagonal 
direction the peak structure for the symmetry broken sys- 
tem approaches the unperturbed form as one approaches 
Ef- On the other hand, scanning along the M^X di- 
rection, the peak structure basically is built up by the 
two bands with intermediate intensity below E^ which 
are shown in Fig. 7b. Therefore one observes a dou- 
ble peak structure at the M-point where the lower peak 
decreases in intensity upon approaching the Fermi level. 
Due to the induced CDW-gap at the FS crossing of the 
homogeneous system k y = (tt, 0.52) there is a shift of the 
spectra to lower energy as can be seen from the right 
panel of Fig (8b). It is also quite interesting to see how 
the CDW-gap occurring at k y — (n, 0.52) evolves when 
one moves along the FS of the unperturbed system. The 
corresponding energy distribution curves are shown in 
Fig. 9. It turns out that the gap first keeps its value for 
small angles up to 8 w 6°. Moving towards the diago- 
nal, the (pseudo)gap moves at energies below the Fermi 
level, so that at the Fermi energy it rapidly decreases and 
vanishes. This occurs at 6 w 20°. In addition it is again 
evident that upon approaching the diagonal direction the 
CDW peak structure evolves into the unperturbed one. 
These findings may correspond to the experimentally ob- 
served shrinking of the FS p9[ 
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FIG. 8. Photoemission spectra for cuts along (a) V < k < 
X, (b) M < k < X. The fermi function has been added as 
a background. Solid lines: Eggbox solution for the same pa- 
rameter set as in Fig. 1. Dashed lines: Homogeneous system. 
The broadening of the f5-functions is 25meV. 

The band structure along M — > X is plotted in Fig. 
7b. Also in this case we observe that at the M-point the 
CDW scattering has induced various bands with inter- 
mediate and low intensity. At the former FS crossing 
(indicated by the arrow) now a gap occurs which is of 
the order of the CDW order parameter \- Although the 
main band (full circles) bends upwards when approach- 
ing the Fermi level there is a FS crossing by the shadow 
band with intermediate intensity. This crossing occurs 
at a slightly larger k y value than in the unperturbed sys- 
tem which mirrors in Fig. 6 as a broadening of the tube 
structure around the M-points. 

In Fig. 8 we show the energy distribution curves along 
the r — > X and M — > X directions. Scanning along the 
diagonal one clearly observes the multiple peak struc- 
ture for the deeply bound states in Fig. 7a. However, 



We now conclude this analysis of the 2-d eggbox CDW 
scattering by considering the relevant issue of the en- 
hancement of the van Hove singularity around the M 
and M points. In particular, as it can be seen in Fig. 5d, 
one finds a generic broadening of the spectra, with disper- 
sionless portions of the spectrum scattered at both higher 
and lower energies with respect to the unperturbed quasi- 
particle band. Moreover, two dispersive shadow bands of 
weak intensity cross the Fermi energy nearly halfway be- 
tween r and M. All these effects contribute to bring 
spectral weight at closer distance from the Fermi level. 
Once the spectral features observed in the present sim- 
plistic mean-field treatment are broadened by the fluctua- 
tive character of the CDW and by the many-body effects 
acting in the real systems, it is not unconceivable that 
our present findings correspond to an extension of the 
dispersionless region of the experimental band structure. 
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FIG. 9. Scan of the photoemission spectra along the un- 
perturbed Fermi surface. The plots are labeled by O which is 
defined in the right panel. Full lines: CDW spectra; Dashed 
fines: Spectra of the unperturbed system. 



C. Disordered eggbox phase 

In the preceeding section we have considered an ef- 
fective one-particle Hamiltonian so that all poles of the 
resulting spectral function correspond to delta functions. 
We now want to extend our results to have some qualita- 
tive idea of a more dynamical description of the problem, 
i.e. the quasiparticle coupling to collective incommensu- 
rate CDW fluctuations. 

As in Ref. [E2] we start from a phenomenological one- 
particle Hamiltonian 



J? = 5^£ fc n Jfc + 5^V(Ri)n i 



(8) 



where the kinetic part is the same as in the preceeding 
section and V(Ri) represents an effective eggbox poten- 
tial given by 



V(Ri)=2V J2 sech (— 



— )sech{— L - 



Un 



i 



) (9) 



The amplitude and broadening of the charge modulation 
is determined by Vq , £ and x n , y n fix the positions of the 
individual 'eggboxes'. In accordance with the ordered 
eggbox modulation we choose a mean charge separation 
of 8a and define x n +i —x n = y n +i — 2/n = 8a + rj where 'a' 
is the lattice constant and r\ is a random number varying 
between —pa < r\ < pa. Taking p — one has again 
an ordered eggbox potential whereas upon increasing p 
(Vmax = 7) long-range charge order is destroyed to an 
increasing degree. The charge amplitude Vq is restricted 
again to values where zero double occupancy at each lat- 
tice site is preserved. The results presented below are 
averages over five random configurations calculated by 
diagonalizing a square lattice with dimension 40x40 in 



real space. Doping is again 6 = 0.2 and the secans-type 
structure is broadened by £ = 2. 

Fig. 10 displays the Fermi surface for a disordered 
eggbox-type charge modulation. The qualitative features 
are the same as in Fig. 6, namely the reduction of in- 
tensity around the M-points. However, the displaced 
shadow bands are no longer visible now. Instead the 
Fermi surface is smeared in k-space i.e. its boundaries 
consist of states with low intensity (< 10%). Around the 
M-points these states completely fill up the tube struc- 
ture. This could correspond to the findings of sequential 
angle-scanning photoemission where diffuse features in 
these regions have been detected J 2(| an d to the suppres- 
sion of the FS observed in ARPES 
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FIG. 10. Fermi surface for a disordered eggbox charge mod- 
ulation. Parameters: S - 0.2, a = -0.45, p - 3, V - 0.12eV. 
The plot is for temperature T = WOK and the energy window 
around Ef has choosen to be 50meV. Intensities: / > 50%: 
full points, 10% < I < 50%: circles, 1% < I < 10%: small 
dots. 

Since we are considering averages over disordered 
charge configurations the energy bands corresponding to 
Fig. 7 now are also smeared along the energy axis there- 
fore removing the fine structure of the CDW gap topol- 
ogy. This can be seen in Fig. 11 where we have plot- 
ted the energy distribution curves along the T — > X and 
M — > X directions in comparison with the homogeneous 
system for two different values of the parameter p which 
defines the degree of suppression of long-range charge or- 
der. This smearing of the individual CDW peaks is of 
course more pronounced for p = 7 in Fig. lib than for 
p = 3 in Fig. 11a. However, scanning along the T — > X 
direction we observe the same feature than for the or- 
dered eggbox modulation, namely the evolution of the 
spectra into a single peak structure upon approaching 
the Fermi level thus recovering the FS segment of the 



10 



quasiparticle. On the contrary, at the M-point the quasi- 
particle peak is now completely suppressed and the spec- 
tra are described by a broad feature extending to very 
low energies. 
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FIG. 11. Photoemission spectra for cuts along V < k < X 
and M < k < X for the same parameter set then in Fig. 
6. The broadening of the (5-functions is 25meV. Dashed lines: 
Homogeneous system. Solid lines: Spectra for the disordered 
eggbox system. The parameter p defing the destruction of 
long-range charge order is p=3 in (a) and p=7 in (b). 



IV. DISCUSSION AND SUMMARY 

In this paper we have investigated the possible con- 
sequences of incommensurate CDW scattering with re- 
gard to the unusual Fermi surface and photoemission fea- 
tures in underdoped bismuth cuprates already mentioned 
above. Our results have been obtained within a scheme, 
which is subjected to some limitations both at the level 
of the starting model and of its treatment. In particu- 
lar, only the physics of CDW modulations has been con- 
sidered, thus neglecting the relevant interplay between 
charge, spin, and Cooper pair fluctuations. Nevertheless 
it is quite interesting and instructive by itself that our 
treatment of a pure CDW symmetry broken system still 



captures some features of the observed spectral proper- 
ties. This indicates that charge-ordering can indeed play 
a role in determining the properties of the cuprates. 

The simplified model considered here has been approx- 
imately treated in a mean-field scheme aiming to capture 
features of well-formed locally ordered stripes in the un- 
derdoped cuprates. In this scheme the fermionic quasi- 
particles do not interact neither among themselves nor 
with the collective CDW fluctuations. Instead the static 
symmetry breaking due to charge modulation produces 
Bragg scattering giving rise to multiple bands and gap 
opening. As a result the bands arising from the mean- 
field description are due to <5-like quasiparticle peaks, 
which obviously do not individually correspond to the 
real experimental features. Nevertheless, it is worth not- 
ing that the multiple Bragg scattering due to incommen- 
suration leads for most fc-points to the appearance of 
multiple quasiparticle peaks [the shadow bands, see Figs. 
(2), (5), (7)]. It is quite natural that the broadening aris- 
ing from the scattering between quasiparticles and the 
scattering between quasiparticles and collective modes 
will mix these peaks thus producing the broad features 
commonly detected in photoemission experiments. Still 
we believe that our analysis provides the location and 
relative position of the spectral features. 

We considered one-dimensional charge modulations 
along the (1,0) and/or (0,1) as well as a two-dimensional 
eggbox-like texture. The one-dimensional solutions could 
account for some enhancement of the van Hove singular- 
ities, but do not seem particularly successful in describ- 
ing the appearance of a pseudogap near the M and M 
points. On the other hand the eggbox case turns out to 
be particularly appealing. The main findings are that 
the 2-d eggbox CDW scattering might account both for 
a flattening of the band dispersion in the r — > M and 
r-tJH directions and for the arising of a leading-edge 
gap around the M, M points leaving finite portions of the 
Fermi surface gapless. As discussed in the introduction, 
this latter quite unusual and non BCS-like behavior of 
the gap has indeed been observed [g9| and suggests that 
the charge fluctuations in the particle-hole channel could 
substantially participate to the spectral features around 
the M, M points. In this regards our finding supports 
the idea that the scattering induced by a charge ordered 
superstructure works in many aspects in the same direc- 
tion than mechanisms involving incoherent pairing cor- 
relations (possibly cooperating with magnetic scattering 
p7[). On the other hand, we also found that a disor- 
dering of the eggbox structures "washes out" the weak 
low-energy features and brings about the disappearance 
of the leading-edge gap due to charge scattering, but 
preserves the feature of Fermi surface only formed by 
disconnected arcs. In this last situation, the pseudogap 
will only appear as a (robust) suppression of the spectral 
weight, which is shifted at higher binding energies. 
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The overall success of the eggbox CDW scattering 
in reproducing some highly non-trivial features of the 
single-particle spectra of the underdoped cuprates sug- 
gests that a dynamical charge pseudo-ordering is a rele- 
vant aspect of the physics of these materials. 
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